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ABSTRACT 


We  describe  a  general  approach  to  the  parametrization  of  linear  combinations  of  Gaussian 
atomic  orbitals,  useful  for  atomic  and  molecular  interactions.  We  use  a  Gaussian  transfcnm 
method  and  Gauss-Legendre  quadratures  to  express  hydrogenic  atomic  orbitals,  with  varying 
effective  charges,  in  terms  of  Gaussian-type  orbitals.  This  procedure  provides  well  defined 
rules  for  calculating  exponent  factors  and  combination  coefficients  of  the  linear  combinations  of 
Gaussians  in  problems  where  nuclear  distances  may  vary  over  large  ranges  during  interactions. 
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1.  INTRODUCTION 


Since  the  suggestion  by  Boys  in  1950,*  the  use  of  Gaussian  functions  as  basis  sets  has 
increased  useful  in  theoretical  chemistry  calculations  of  electronic  structure.  Gaussian  functions 
have  been  found  very  useful  in  the  efficient  calculation  of  multicenter  electronic  integrals. 
Gaussian  basis  sets  have  been  parametrized  to  optimize  the  calculation  of  atomic  and  molecular 
energies,  in  work  that  has  continued  until  voy  lecently.^*^  Gaussian  functions  also  arise  when 
one  introduces  the  Gaussian  transform  of  atomic  orbitals  (AOs).  In  this  case,  the  transfonn 
integral  is  eventually  discretized,  to  obtain  a  set  of  parametrized  basis  functions.  The  choice 
of  discretization  has  again  been  guided  by  interest  in  the  calculation  of  energies  of  stationary 
states.^’*^*** 

Basis  sets  are  also  needed  in  studies  of  atomic  and  molecular  interactions,  where  the  main 
concern  is  the  calculation  of  probabilities  of  electronic  transitions,  exdtadtm  cross  sections,  and 
thermal  rates.  We  have  used  atomic  bases  in  recent  studies,  with  a  time-dependent  description 
which  construct  molecular  orbitals  as  combinations  of  AOs  with  time-dependent  coefficients.*^ 

In  the  case  of  atomic  or  molecular  interacticms,  where  distances  and  bond  angles  change 
over  largp  values,  it  is  important  to  work  with  AOs  having  the  correct  asymptotic  behavior.  This 
is  given  by  radial  functions  which  decrease  exponentially  with  the  distance  from  an  electron 
to  an  atomic  cote,  at  a  rate  depending  on  the  effective  charge  of  the  core.  This  behaviour  is 
ccnrectly  given  by  hydrogenic  atomic  orbitals  (HAO’s).  Their  orbital  exponents  must  change 
to  reflect  varying  effective  nuclear  charges  as  nuclei  or  cores  move  and  exponents  change  with 
time.  The  aim  of  this  contribution  is  to  provide  a  general  way  of  expressing  these  HAO*s  in 
terms  of  Gaussian  atomic  orbitals  (GAO’s),  with  exponents  and  coefficients  directly  determined 
from  the  varying  effective  charges.  Our  criterion  for  choosing  the  Gaussian  basis  sets  is  that  they 
should  properly  describe  the  AOs  at  distances  from  atomic  cores  that  are  important  in  dynamical 
phenomena,  instead  of  emphasizing  the  energetics  of  stationary  stares. 

The  angular  part  of  the  basis  fimction  can  be  represented  in  terms  of  Qtherical  harmonics 
or  in  terms  of  integer  powers  of  the  Cartesian  cocmlinares.  In  this  case  a  primitive 
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ftnctkm  may  be  wrinra  u 

*»(*» »» *)  *  »"•  e*P(-a*^)»  (1) 

and  a  contracted  Gaussian  basis  funcdtm  is  represented  by’ 

x(*»  y.  *)  *  **'  y"*  **'  Ci  e*p(-«iT*).  (2) 

• 

We  have  recently  used  the  tune-dependent  SchrOdinger  equation  and  a  related  time-depuident 
variational  principle,  to  obtain  complex  time-dependent  hydrogenic  atomic  orbitals  (CTDHAO’s) 
resulting  from  their  interacticm  with  a  moving  positive  charge;^^  this  procedure  leads  to  optimized 
trial  AOs  with  time-dependent  orbital  expmwnts.  Here  we  use  the  Gaussian  transform,  which  also 
rq^)ears  in  a  generator  coordinate  crmtext.*’  to  express  these  CTDHAO’s  in  terms  of  Gaussian 
functitms.  This  allows  us  to  construct  time-depemfent  Gaussian  basis  sets  at  any  given  time. 
The  CTDHAO’s  have  the  form 

where  the  time-dependence  appears  through  the  complex-valued  parameters 

Cn/(<)  *  Kii*)  +  W 

We  have  found  from  the  calculaticm  of  the  time  evolution  of  the  complex  parameters'^  that 
K  <  A  for  charges  moving  at  the  thermal  and  hypertheimal  velocities  of  interest  in  chemical 
qrplications;  therefore  we  use  the  approximation 

m  «  m  (5) 

for  die  mqKmential  parameters.  From  now  cm  deal  with  HAO’s  whose  varying  effective  charges 
have  real  values.  Ws  introduce  a  convenient  change  of  integraticHi  variable  in  their  Gaussian 
transform  integral,  befrae  discretizing  the  variaUe  by  means  of  the  Gauss-Legendre  quadrature.*’ 
^  can  control  the  accuracy  of  our  Gaussian  fits  by  increasing  the  number  of  quadrature  points 
in  a  systenoatic  way.  This  allows  us  to  obtain  excellent  fits  to  the  HAOs  without  requiring  least 
square  fits  or  energy  optimizathm. 
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This  ]My)er  is  ocguiz^  as  follows.  In  sections  2  and  3  we  show  how  to  express  the  HAO*s  in 
tenns  of  a  Gaussian  basis  set  of  given  size.  In  section  4  we  give  expressions  for  onhogooalizing 
die  Ofbitals  by  the  Oram-Schmidt  inocedure,  and  in  sectum  S  we  express  them  in  tenns  of 
cartesian  fl^nfi****- 

2.  HYDROOENIC  ATOMIC  ORBITALS 

Hie  ndial  pans  of  the  Is.  2s.  and  2p  HAO's  ate 

Rl.(r:A,.).2A,.*e-*".  (6) 

Aj.)  =  •^A  J(1  -  ^)e- V,  (7) 

and 

Rsalr;  As,)  -  ^A„*  re-¥,  (8) 

where  the  expcmential  parameten  may  change  over  time  due  to  nuclear  motions. 

3.  GAUSSIAN  BASIS  SETS 

We  next  express  the  HAO’s  in  terms  of  Gaussian  transforms.  This  procedure  allows  us  to 
choose  the  best  set  of  Gaussian  functions  for  example  at  any  given  time  during  the  coUisicm 
process. 

A.  The  Is  orbital 

hi  the  Is  case.  Eq.  (6)  becomes 

Rl.(ri(l.)»2A„*.-»-'.  (9) 

In  order  to  eliminate  A^i  in  the  pre-exponential  term  we  make  a  rescaling;^ 

Pnt  *  A*|  r.  (10) 

Now  the  normalized  radial  part  of  the  orbital  is 

(11) 
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which 


<lfi.|5i.>  -  - 1.  (12) 

it 

The  Gausdan  tunsform  of  1Ea«(pi«) 

oo 

0 

where  Ac  weight  function  7i.(») » 

7,.(Zf)  =  (<r5*r»  «-*,  (M) 

and  the  geirerator  function  G(Zk,pia)  is 

G(s.n.)  =  03) 

wiA  a  being  the  generator  coordinate,  i  e..  the  Gaussian  exponential  parameter. 

We  are  interested  in  numerically  evalutuing  Eq.  (13)  wiA  a  small  number  of  Gaussian 
functions.  So  we  malre  Ae  following  change  of  variables 

V  * 

which  transforms  the  mtegraticm  mterval  fmn  (0,oo)  to  (0,1). 

Eq.  (13)  then  becomes 


Using  the  Gauss-Legendre  quadrature^*  Ae  mtegral  m  (17)  becomes  a  summation  and  we 


where  die  weights  toj  and  die  abscissas  q,  give  the  set  of  exponent  factors  and  combination 
coefficients.  Ws  have  used  Gaussian  functions  duou^iout  this  paper. 


Eq.  (18)  has  two  singular  values  at  7"0  and  i7«l.  In  the  present  case,  where  we  are  using 
a  small  number  of  points,  this  is  not  a  problem  because  the  limit  points  are  reastmably  far  from 
0  and  1,  reqiectivdy. 
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B.  The  2s  orbital 

For  die  2s  cubital,  Eq.  (7),  we  have  after  the  rescaling  of  Eq.  (10)  and  nonnalizadtm  that 


K2.(«.)  =  ^(1  -  ^)e-% 

(19) 

which  has  the  fdlowing  Gaussian  transform^ 

00 

0 

(20) 

where 

7i.(S)  =  (8ir5‘r*(2-^)e‘*. 

(21) 

and 

G(3,(>J.)>=  e-*^. 

(22) 

Applying  the  change  of  variables  of  Eq.  (16)  to  Eq.  (20)  we  have 

1 


(-Ini7)i  '^*''l61ni7 


which  becomes,  after  application  of  the  Gauss-Legendxe  quadrature 

If  tn  ,  n\ _ \  _  2 


v  /  \  ^  (2  +  21niyi)  ,  P2i  , 


C.  The  2p  orbitals 

Fbr  the  2p  orbitals,  Eq.  (8),  the  equivalent  of  Eq.  (19)  is 


The  Gaussian  transform  is*^ 

eo 

*=  J 7ij»(o)G'(q,  p2f)dSf 
0 


(23) 


(24) 


(25) 


(26) 


s 


'V.  '  ’  ■ 


^uhat 

and 


7j,(S)  -  (S4irB*)-J  (i  -  1)  e-A, 

G(®i«j)“  e*®^. 


Mth  die  chan^  (tf  variables  in  Eq.  (16),  Eq.  (26)  becomes 


or 


(-l~21niyi) 

(-In^O^ 


exp( 


P2r 


IGlnijt 


)dr). 


...  (-l-21niyi) 

^  (-InnOi 


exp( 


P2,^ 

ISlnt;,- 


). 


(27) 

(28) 

(29) 

(30) 


Widi  the  use  of  Eq.  (10)  we  can  write  Eqs.  (18),  (24),  and  (30)  as 

s 

=  (31) 

ml 

where  the  linear  coefficients  c^ij  and  the  exponents  tVi  (p^en  in  Thble  I)  are  expressed  in  terms 
of  the  weights  and  abscissas  of  Gauss-Legendre  quadrature,  which  are  given  in  Table  IL 
Rni  (>*;  Ant),  Eqs*  (6)-(8),  can  be  obtained  from  ^  (pnt)  by  noting  that 


Rut  (**i  Ant)  .  J 
Snt  (Pni)  ‘  * 


(32) 


4.  CARTESIAN  GAUSSIAN  REPRESENTATION 

Because  they  have  been  preferred  in  molecular  calculations,  we  are  going  to  express  our 
sphorical  Gaussian-type  oibitals  in  terms  of  cartesian  Gaussian  functions. 

We  start  out  re^xpressing  the  HAO’s,  Eq.  (3),  as 

0«im(r,  t)  *  /2ni [r;  Ani(t)]r~V* Jlm(^,  v),  (33) 

where  the  solid  spherical  harmonics.  r'}^m(^>  V’)*  correspond  to  linear  combinations  of  the  integer 
powers  of  the  cartel  coordinaies. 
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In  the  case  of  s  oibitals.  M)  and  the  cartesian  and  spherical  representations  are  the  same. 
In  the  case  of  2p  otitals  we  have 

w^iere  k  stands  for  x,  and  x. 

Applying  the  rescaling  of  Eq.  (10)  to  the  "radial**  part.  Eq.  (34)  becomes 

^Pk  (**.  0  - 


which  has  the  fcdlowing  Gaussian  transfonn 


OD  2 

0 

Applying  the  change  of  variables  (16)  and  the  Gauss-Legendre  quadrature  to  Eq.  (36)  we 
have 


V»2p»(r,<) 


1 


(37) 


If  we  express  Eq.  (37)  like  we  did  in  Eq.  (31)  we  have  that 

V^P*(r.  ^  cxp(-a2,,,T*) .  (38) 

tael 

where  «nd  axpj  are  given  in  tidtle  I. 

5.  GRAM-SCHMIDT  ORTHOGONALIZATION 

The  exact  hydxogen>like  orbitals  are  onhononnal,  L  e.» 

{^mlm  I  (39) 

However  when  we  calculate  Eq.  (39)  using  six  Gaussian  functions  we  obtained  0.015976, 1.0015, 
0.892141,  and  0.942533  for  die  ls>2s,  Is-ls,  2s-2s,  and  2p-2p  overliqis,  respectively.  Therefore 
we  decided  to  orthogonalize  tte  orbitals  with  the  (}ram-Schmidt  inocedure.*’ 
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(40) 


yfh  begin  with  die  nonnalization  of  ICi«.  If  we  take 


and  impose  that 


then  we  have  that 


(pi«)  ®  NijRi,(pi,), 


(ErPi."’) = /  ipi.  rj'.cpu)’ = 1, 


where 


Ni.= 


'(Ru^Ru) 


K  M 


(Ei.|E..>  =  -i.EE7^ - 7^ 

(43) 

r— J_4.— 

Next  we  want  to  orthogonalize  the  2s  oifaital.  To  do  this  we  take  the  new  orthogonal 
unnoimalized  2s  orbital,  '^g{p2t)i  ^ 

Kifu)  =  Sj.(/>2.)  +  ^u\pu).  (44) 


The  lequiiviiieat  that  be  oithogonal  to  Ru\pi*)  allows  us  to  detennine 

^=-(Rs.ii5r>= 


-N  1  w.  Mii(2  +  21nUi),  1  ,  1  ,-t 

**  2\/2» " "  4(-lni^,)  16(-lni7,)^ 

In  the  present  case  ^«-0.015963. 

Now  we  normalize  the  2s  orbital  by  expressing  it  as 


^2}  \p2»)  *=  N2,1^(P2«)- 


Normalizing  '  to  unity  gives  us 


N,.= 


(Ra^l^) 
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where 


N  N 


«  4'W 


tt>,(2  -h  2 In )?,)  Wj(2  -t-  2 In rij) ,  1 


yf-^j  16(-lni?<)  16(-lni?j) 


Since  we  cannot  represent  ^«(p2«)«  Eq.  (44),  as  a  single  summation,  we  have  instead 
shown,  in  Ihble  m,  the  normalized  ^2«(/>2«) 


^\n.)  =  N,. !,.(«.), 


whoe 


N2,= 


^<R2.1R2.) 


In  the  case  of  2p  orbitals  we  just  nonnalize  them  to  unity  because  their  spherical  harmonics 
make  them  orthogonal  to  s-orbitals.  So  the  nonnalized  radial  part  is 


\p2p)  *  N2pS2p(/>2p)» 


where  (by  the  normalization  condition) 


'<R2p|R2p> 


(^2p(^2p)  = 


if  N 


_  ^  V'  ^  tgi(-l-21nT?.)u>,(-l~21ni?,)f  1  ^  1  (53) 

v'-iniy.  v^-ln7y  46(-lni7i)  16(-lni/>)' 


In  the  case  or  cartesian  2p  orbitals  we  normalize  the  whole  orbital  to  unity ,  i.  e., 

»  N„.  g_^exp(^) , 

wtere  (by  the  normalization  condition) 
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integrals  in  molecular  calculations.  An  altrnnadve  is  to  orthonormalize  the  oibitals.  The  Gram- 
Schmidt  procedure  was  found  to  be  ctmvenient  for  this  purpose;  it  can  be  implemented  with 
the  expressions  we  presented  in  Section  S,  and  require  only  a  knowledge  of  the  Gauss-Legendre 
quadrature  points.  The  small  calculated  value  of  the  ^  coefficient  there  suggests  that  forcing 
more  strict  orthogonality  will  have  only  a  small  effect  on  molecular  densities. 

The  figures  show  that  the  behavior  of  the  physical  HAO’s  can  be  well  represented  over 
wide  distances,  as  required  for  molecular  interaction  studies.  The  largest  discrepancy  is  found 
near  the  origen  for  the  cusp  behaviour,  as  expected.  This  region  of  the  atomic  distances  is 
however  inaccessible  for  thermal  and  hypeithermal  collision  velocities  so  that  an  error  in  that 
region  would  not  have  a  noticeable  effect  on  collisional  properties  such  as  cross  sections.  The 
long  range  behaviour  of  the  HAOs  is  very  well  reproduced  and  indicates  that  the  constructed 
orbitals  would  have  the  correct  interatomic  overlaps  required  for  the  description  of  electronic 
rearrangement  in  molecular  interactions. 


One  of  the  authors  (HFMDC)  would  like  to  acknowledge  D.  BekSi6  and  K.  Runge  for 
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Table  I 

Linear  coefficients  c,-  and  exponents  a,-  for  the 
Gaussian  expansion  of  the  HAO’s. 


Orbital 

Ci 

ai 

Is 

7?(-b;)T 

\  * 

,^|JU 

41a 

2s 

2p 

1  * 

^29 

16  la  rii 

1  ^ 

\  2 

ISlaiji 

4 


i 


Table  H 

Weights  Wi  and  abscissas  rji  for  the  Gauss- 
Legendre  quadrature  with  N=6. 


Wi 

»?• 

0.085662 

0.033765 

0.180381 

0.169395 

0.233957 

0.380690 

0.233957 

0.619310 

0.180381 

0.830605 

0.085662 

0.966235 
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f 


Table  HI 

Numerical  values  of  linear  coefficients  Ci 
and  exponents  a,-  for  normalized  orbitals  with 

A/i/— 1.0.. 


i 

ESBIHli 

C2p,i 

(carl) 

*^P,i 

1 

0.052469 

0.073783 

-0.093888 

0.018446 

0.063778 

0.002621 

2 

0.152629 

0.140804 

-0.406076 

0.035201 

0.081931 

0.007626 

3 

0.268416 

0.258861 

0.006884 

0.064715 

0.052614 

0.013412 

4 

0.381074 

0.521758 

0.148709 

0.130439 

-0.003344 

0.019041 

5 

0.472074 

1.346976 

0.288046 

0.336744 

-0.062462 

0.023588 

6 

0.521129 

7.278401 

0.377034 

1.819600 

-0.102124 

0.026039 
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Figure  Captions 


Figure  1  —  Comparison  between  the  exact  radial  part  of  the  Is  normalized  hydtogenic 
orbital  and  the  corresponding  numerical  one  obtained  in  this  woik  with  6  Gaussians. 

Figure  2  —  The  same  as  figure  1  for  the  2s  oibitaL 
Figure  3  —  The  same  as  figure  1  for  tl^  2p  orbital. 

Figure  4  —  Comparison  between  the  exact  “radial*’  part  of  the  cartesian  2p  orbital,  given 
by  (3/4x)  -1/2  ^-1  corresponding  numerical  values  obtained  in  this  work  with  6 

Gaussians. 
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